Evidence of electro- active excitation of the spin cycloid in TbMnOs 



o ■ 
(N : 



I 

O 

o 



> 
oo 



o 

O 



X 



A. M. Shuvaev/ V. D. Travkin,^ V. Yu. Ivanov,^ A. A. Mukhiii,^ and A. Pimenov^ 

^ Experimentelle Physik 4, Umversitdt Wiirzburg, D-97074 Wiirzburg, Germany 
^General Physics Institute, Russian Academy of Science, 119991 Moscow, Russia 

(Dated: January 7, 2010) 

Terahertz electromagnetic excitations in the multiferroic TbMnOa at the field-induced magnetic 
transition are investigated for different orientations of the magnetic cycloid. In addition to the elec- 
tromagnon along the a-axis, the detailed polarization analysis of the experimental spectra suggests 
the existence of an electro-active excitation for ac electric fields along the crystallographic c-axis. 
This excitation is possibly the electro-active eigenmode of the spin cycloid in TbMnOs, which has 
been predicted within the inverse Dzyaloshinskii-Moriya mechanism of magnetoelectric coupling. 



tion must be parallel to the c-axis (f ||c) and the electro- 
magnons should be excited for ac electric fields parallel to 
the a-axis (e||a). Accordingly, for the ab-plane oriented 
cycloid the polarization P\\a and electromagnons for e|jc 
should be observed. This basically corresponds to the in- 
terchange of the a and c-axes in the excitation conditions. 
Although the orientation of the static electric polariza- 
tion has been confirmed in many experiments J5i, l24i . 1251] , 
the excitation conditions for the electromagnon contra- 
dicted the predictions of the model: the selection rules 
for the excitations remained the same (i.e. Sjl a) inde- 
pendently of the orientation of the cycloid llHlOj. In 
order to account for this contradiction, a model based 
on Heisenberg exchan ge c oupling between spins has been 
proposed recently [H, 3, i^]. According to this model, 
the structural peculiarities of orthorhombic multiferroic 
manganites lead to the selection rules e\\a regardless of 
the orientation of the spin cycloid. 



Multiferroics represent an intriguing class of materials 
in which magnetism and electricity are strongly coupled. 
This magnetoelectric coupling leads to a mutual influ- 
ence of magnetic and electric ordering, which results in a 
rich new physics and in various possible applications [JJ- 
0]. Recently, a new class of multiferroics has attracted 
enormous interest because the ferroelectric polarization 
in these systems is induced by a cycloidal ordering of the 
magnetic moments [Hl-Q- The magnetic and electric or- 
der in these materials are accompanied by coupled spin- 
lattice vibrations ^] showing strong electric dipole activ- 
ity and termed electromagnons 0, [l3] ■ Initially detected 
in GdMnOs and TbMnOa ll'], the electromagnons ap- 
pear to exist in cycloidal magnetic phases of many differ- 
ent multiferroics [l^ - [l9l |. In most cases, distinct struc- 
ture of these excitations is seen including two or more 
separate modes. For example, in TbMnOa one can ob- 
serve at least 3 excitations: two with the energy around 
~ 2 meV (10-20 cm~^) and the third one at ~ 8 meV 
(60 cm~^). 

Ferroelectric polarization in spiral magnets has been 
successfully explained on the basis of the inverse 
Dzyaloshinskii-Moriro (DM) coupling between the mag- the nature of the lower energy elecromagnons close to 
netic moments ji^-lli]. Within this model a spiral spin 20 cm'^ remains uncertain. It seems reasonable to re- 
call the statements of the DM model and to assume that 
the low-frequency electromagnons are somehow related 
to the excitations of the s pin cyc loid observed by inelastic 
neutron scattering (INS) [29|-|31 [because the characteris- 
tic frequencies coincide closely 29, 32]. This assumption 
is supported by a recent observation of the magnetic exci- 
tation channel for electromagnons [55] which can be also 
seen as antiferromagnetic resonances (AFMR) within cer- 
tain excitation conditions. However, in order to prove 
this assumption, the predicted excitation conditions for 
the eigenmodes of the cycloid should be found experi- 
mentally. Most specifically, one should expect the ex- 
citation conditions along the c-axis if the spin cycloid is 
oriented in the a6-plane. Such excitation conditions were 
not observed up to now |2^, [s^ . A possible reason for 
this fact is the weakness of the dielectric contribution of 
the spin modes within the DM mechanism. In order to 
resolve this experimental difficulty, TbMnOs seems to be 
an ideal candidate, because the magnetic cycloid can be 



The Heisenberg exchange model assigns the high en- 
cm"^ and 
but 



ergy electromagnons observed between 60 
80 cm~^ to the zone edge magnons [l 



26-2 



structure lowers the symmetry of the system leading to 
an effective magnetoelectric coupling term and an elec- 
trical polarization P in the form 



X Si X S j 



(1) 



Here Si and Sj are the neighboring spins and is the 
lattice vector connecting them. Consequently, it was rea- 
sonable to assume the same mechanism to explain the 
existence of the electromagnons 23|. In this model the 



electromagnons are the eigenmodes of the spin cycloid 
which can be excited by electric component of the elec- 
tromagnetic wave However, after initial success of 
this explanation a number of experimental results con- 
tradicted the predictions of the model. The basic exam- 
ple can be given by TbMnOa at low temperatures: the 
DM model predicts that the excitation conditions for the 
electromagnon should be tied to the magnetic cycloid, 
i.e. for the bc-plane oriented cycloid the electric polariza- 
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rotated between a6-plane and &c-plane in external mag- 
netic fields. Thus the ultimate experiment must follow 
the c-axis response at the magnetically induced rotation 
of the cycloid. 

In this work we have carried out detailed polarization 
analysis of terahertz excitations in TbMnOs. Special at- 
tention has been payed to the selection rules of different 
excitations upon the rotation of the spin cycloid from be- 
ta a6-plane in external magnetic field _B||5-axis. It was 
found that a new excitation arises in the high-field phase 
with excitation conditions e||c and h\\a. We argue that 
this excitation could not be explained by purely magnetic 
contribution but carries substantial electric component. 

The transmittance experiments at terahertz frequen- 
cies (3 cni^^ < ly < 30 cm~^) have been carried out 
in a Mach-Zehnder interferometer arrangement 33, |33| 



which allows measurements of amplitude and phase shift 
in a geometry with controlled polarization of radia- 
tion. The experiments in external magnetic fields up 
to 8 T have been performed in a superconducting split- 
coil magnet with polypropylene windows. Single crys- 
tals of TbMnOa have been grown using the floating-zone 
method with radiation heating. The samples were char- 
acterized using X-ray, magnetic, dielectric and optical 
measurements 0, [3^1 . The results of these experiments 
including the magnetic phase diagram are closely similar 
to the published results [j^ ]. 

The paramagnetic phase in TbMnOa above 40 K is 
followed by a sinusoidally-modulated antiferromagnetic 
state which transforms below 28 K into the cycloidal spin 
structure oriented within the &c-plane [l^ [stL (inset in 
Fig. [1]). According to the symmetry analysis 
Eq. ([T]), static electric polarization along the c-axis is 
allowed in the low-temperature phase. In external mag- 
netic fields t he sp in-cycloid rotates from &c-plane towards 



the a&-plane [25|,|3H. Correspondingly, the electric polar- 
ization rotates from P|| c-axis to the P||a-axis [sl. [jj. [25| . 

Dynamic experiments on TbMnOs in external mag- 
netic fields reveal that the c-axis properties are indeed 
sensitive to the orientation of the cycloid. The exam- 
ples of such changes are shown in Fig. [1] which represents 
the terahertz properties of TbMnOa in external magnetic 
field B\\b inducing the transition from 6c-plane to ab- 
plane oriented cycloid. The magnetic field scans in these 
experiments were made in the geometry with e||c and 
h\\a and at T = 10 K. The data are represented as re- 
fractive index n+iK — ^feji as both electric and magnetic 
contributions could be mixed in this experimental geom- 
etry. The transition from fee-plane to o6-plane cycloid in 
the high magnetic fields is clearly seen. Already at this 
point it is evident that the observed changes are strongly 
frequency dependent. Here the data at 4.7 cm~^ is in- 
fluenced by a Tb-mode around 5 cm~^ [s^] which disap- 
pears in the high-field phase with the a6-plane cycloid. 
This leads to a substantial decrease of the absorption 
(k(4.7 cm^^)) and reveals a bit complicated structure in 
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FIG. 1: Dependence of refractive index (n = Re{y/e]T), up- 
per panels) and absorption coefficient (k — Im{y/ejl), lower 
panels) in TbMnOs on magnetic field _B||6-axis at various fre- 
quencies on crossing the phase transition from the be- to the 
a6-oriented magnetic cycloid. Polarization of incident wave is 
e||c, h\\a, where e and h are electric and magnetic ac fields 
of the electromagnetic wave. The inset shows B-T phase dia- 
gram of TbMnOs for B\\b "ab" and "6c" denote a6-plane 
and 6c-plane oriented cycloids, respectively, PM - paramag- 
netic and IC - sinusoidal phases. 



refractive index below 10 cni^^. The changes observed at 
4.7 cm~^ can be well understood assuming a suppression 
of a Lorentzian mode situated between 5 and 6 cm~^. 

Three higher frequency scans (15.8, 20, and 28 cm^^) 
in Fig. [T] show more systematics, and can be reduced 
22| and to a growth of the absorption mode around 20 cm~^ in 



the phase with the a6-plane cycloid. This is a typical 
behavior for a Lorentz oscillator which appears close to 
20 cm~^ simultaneously with the a&-plane cycloid. In- 
deed, strong additional absorption arises near the fre- 
quency ~ 20 cm~^ and is substantially reduced above 
and below this frequency (lower panels, 28 cm"-'^ and 
16 cm^^, respectively). In addition to a strongly in- 
creased absorption around the resonant frequency, the 
refractive index in Lorentzian model may in rough ap- 
proximation be simplified to the following statements: 
n(j/) w rioc + An for v < i/q and n{u) « Uoo — ^niyt^jv)^ 
for V > vq. Here Uaa is the contribution from the 
high-frequency processes. An refiects the strength of the 
Lorentzian mode, and k, 20 cm~^ is the resonance 
frequency. Therefore, after the growth of the Lorentzian 
(An = — > An ^ 0) at the magnetic transition we 
expect an increase of the refractive index below the 
resonance frequency and a decrease above it. This is 
roughly the behavior of the refractive index, observed 
at 15.8 cm~^ {v < i^q) and 28 cm~^ (t/ > i/q). The 
main changes, seen in Fig. [l] appears at the magnetic 
transition and are therefore related to the rotation of 
the magnetic cycloid. However, better understanding of 
the underlying processes can be obtained within a direct 
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FIG. 2: Transmittance spectra of TbMnOa in external mag- 
netic fields B\\b for different experimental geometries. Sym- 
bols are experimental data and solid lines are fits with Lorentz 
oscillators as discussed in the text. 



analysis of the spectra in the relevant frequency range as 
presented below. 

Figure [5] shows the field dependent spectra for two dif- 
ferent geometries of the experiment. The thicknesses of 
the samples are similar for both orientations: 1.24 mm 
(upper panel) and 1.33 mm (lower panel), respectively. 
The spectra in the lower panel with e||c and h\\a cor- 
respond well to the known results [l^, [s^] and show a 
mode at about 21 cm~^ which appears after the field in- 
duced reorientation of spin cycloid to the a6-plane. Based 
on the weakness of this mode, both in Ref. 26] and in 
Ref. [32] it has been concluded that the mode around 
21 cm^^ is of purely magnetic origin and represent an 
antiferromagnetic resonance of the magnetic cycloid for 
A||a-axis. Indeed, the strength of this mode (Ae ^ 0.05, 
assuming electric origin) is extremely weak cornpared to 
electromagnon observed for e\\a (Ae ~ 2) [^, hj. The 
mode in Fig. [2] is observed for the a6-plane cycloid and 
within h\\a excitation conditions. Tracing this mode back 
into the 6c-oriented cycloid in zero external fields, it can 
be expected to originate from the excitation conditions 
h\\c. (This corresponds to the interchanging of the a-axis 
and c-axis). Indeed, an AFMR mode excited for h\\c of 
the similar strength has been observed around 21 cm^^ 
in the phase with the 6c-plane cycloid [32l |. 

A careful comparison of both panels in Fig. [3] reveals 
interesting difference between two excitation conditions. 
The strength of the mode in the geometry where e\\b is 
roughly the half of that where e|| c. This strongly suggests 
that for geometry in which e||c the electric dipole contri- 
bution is indeed measurable and represent the previously 
unobserved e||c counterpart of the electromagnon. These 



g 10 



. Fi 


a.e 




B.OT 


. Fi 


a.e 




B.5T 


• S 


a,e 




B.7T 


t B 


a,e 


b 


B.8T 




X (cm ) 



FIG. 3: Linear absorption coefficient a = ATrK\~ in TbMnOa 
for different polarizations showing the emergence of an addi- 
tional mode in the phase with the afe-cycloid (green triangles 
and blue crosses). Symbols are experimental data calculated 
from transmittance spectra as described in the text. Lines are 
model calculations using Lorentz oscillators with the same 
parameters as in Fig. [51 The data are shifted vertically to 
eliminate background absorption in different samples. 



results agree well with the original explanation of the 
electromagnons as electrically active eigenmodes of the 
cycloidal structure [9,, ^3] . 

In order to make the discussion quantitative, the ex- 
perimental spectra in the upper panel of Fig.[5]were fitted 
with magnetic Lorentz oscillator. If we now take the pa- 
rameters of the mode from the geometry with e\\b and 
plot the expected transmittance spectra for the geome- 
try e|jc we obtain the absorption value which is too weak 
compared to the experiment (the only" curve in the 
lower panel of Fig. The only possible explanation is 
that this mode has distinct non-zero electric contribution 
along the c-axis. The actual fit for this geometry was 
obtained by taking parameters of the magnetic oscillator 
from e\\b, h\\a geometry and adding an electric oscillator 
with the same resonance frequency — 20.7 cm^^ and 
line width 7 — 4.9 cm~^ as the magnetic one. The rea- 
soning behind this assumption is that both contributions 
are electric and magnetic parts of the same eigenmode 
of the spin cycloid. The strengths of both components is 
given by Afi = 0.0038 and Ae = 0.05, respectively. 

Figure [3] shows the linear absorption coefficients a = 
47rKA~^ which allow direct comparison of different ex- 
perimental geometries. Symbols were calculated from 
transmittance spectra T using the expression T — (I ~ 
i?)^ exp(— ad), where R = \{y/FfI — l)/{^/FfI+ l)p is the 
refiectance on the boundary between the air and the sam- 
ple, and d is the sample thickness. Solid lines are model 
calculations using the same parameters as in Fig. [2j One 
can see again that the mode at 21 cm~^ is stronger in the 
case of e||c and h\\a excitation. (A weak narrow mode 
seen close to 22 cm~^ for e||5, h\\a geometry is possibly 
due to impurities in the sample. The strength of this 
mode is at least an order of magnitude smaller than the 
strength of the broad mode and doesn't change the over- 
all picture.) 
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The mode intensity for the "main" e\\a electromagnon 
(Aeq ~ 2 [ll[) is about 40 times stronger than electric 
contribution along the c-axis (Ae^ — 0.05) observed in 
the present experiment. The large electromagnon ab- 
sorption along the a-axis was one of the challenging ques- 
tions in explaining its origin. The relatively weak static 
electric polarization doesn't fit well with the large dielec- 
tric absorption of electromagnon if both are caused by 
Dzyaloshinskii-Moriya interaction |26|. On the contrary, 
the Heisenberg exchange mechanism (isL 26 . 27| seems to 
explain well the intensities of at least the high-frequency 
electromagnons above 50 cm^^. In this model the edge- 
zone magnon couples to alternating orthorhombic distor- 
tions at oxygen sites via symmetric Heisenberg exchange 
interaction. This leads to the coupling of the zone edge 
magnon to homogeneous electric fields along the a-axis. 
As the symmetric interaction is much stronger than the 
relativistic DM coupling, the hybridized electromagnon 
has enough strength to explain the experimental inten- 
sities for e\\a. Much weaker DM component can- 
not be seen in this experimental geometry because of the 
dominance of the intensity induced by the Heisenberg ex- 
change coupling. On the contrary, rotating the magnetic 
cycloid towards the a6-plane, both contributions can be 
well separated experimentally. The Heisenberg exchange 
part remain oriented along the a- axis, as confirmed by 
different experimental groups EEiliilll. The weak 
DM electromagnon rotates with the cycloid and can be 
clearly observed in the present experiment as electric con- 
tribution along the c-axis. 

One remaining question is: why we observe only one 
mode in the high-field phase? The probable reason is 
that one of two modes is too weak and is not seen in 
the spectra. This argument is supported by recent in- 
elastic neutron scattering experiments (3ll |. In these ex- 
periments the modes of the a6-plane spin cycloid have 
been investigated. Although this a&-plane orientation 
has been achieved using external magnetic field along 
the a-axis, the comparison to the present results is still 
very instructive. It has been observed that the excita- 
tions of the a6-cycloid are dominated by a strong mode 
at 2.25 meV [sij. This frequency corresponds well to the 
excitation at 21 cm~^, seen in Figs. 12131 

In conclusion, we performed detailed polarization anal- 
ysis of the electric and magnetic excitations in TbMnOa 
in the high-field phase where spin cycloid rotates from 
he- to a6-plane. The observed excitation at 21 cm~^ 
could not be described by purely magnetic contribution 
as was suggested previously. We argue that this excita- 
tion is the missing electro-active eigenmode of the spin 
cycloid. The weakness of this mode is in agreement with 
the Dzyaloshinskii-Moriya contribution to the dynamical 
magnetoelectric coupling in TbMnOa , which represents a 
relativistic counterpart to the Heisenberg exchange mech- 
anism. 
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